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Abstract. Complex light metal hydrides possess many properties which make them 
attractive as a storage medium for hydrogen, but typically, catalysts are required to 
lower the hydrogen desorption temperature and to facilitate hydrogen uptake in the 
form of a reversible reaction. The overwhelming focus in the search for catalysing 
agents has been on compounds containing titanium, but the precise mechanism of 
their actions remains somewhat obscure. A recent experiment has now shown that 
fullcrcnes (Cgo) can also act as catalyst for both hydrogen uptake and release in lithium 
borohydride (LiBH^). In an effort to understand the involved mechanism, we have 
employed density functional theory to carry out a detailed study of the interaction 
between this complex metal hydride and the carbon nanomaterial. Considering a 
stepwise reduction of the hydrogen content in LiBH4, we find that the presence of Cqq 
can lead to a substantial reduction of the involved H-removal energies. This effect is 
explained as a consequence of the interaction between the BH~ complex and the Ceo 
entity 
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1. Introduction 

Hydrogen possesses many features which make it a highly attractive option as a future 
energy carrier [TJ |2[ El 0] . While such a transition is certainly not going to be easy, it is 
nonetheless inevitable to explore feasible alternatives to our fossil fuel resources which 
are finite and might even become exhausted within this century, and which furthermore 
involve carbon dioxide emissions linked to a harmful global climate change. Before one 
could even begin to realize such an envisioned hydrogen economy, a number of formidable 
technological problems need to be solved first. One of these issues is hydrogen storage 
[H El El H] . Pressurized gas tanks and cryostatic storage of hydrogen as a liquid do not 
seem ideal for mobile (vehicular) applications. Rather, storage of hydrogen inside solid 
materials appears to be the most feasible solution. 

In this context, complex light metal hydrides are very promising as a suitable 
hydrogen storage medium. On the plus side, the gravimetric hydrogen content in 
these materials can be quite large, but a significant drawback is the generally high 
temperature required for hydrogen desorption and the fact that the reactions are often 
not reversible, i.e., hydrogen can be released, but the direct rehydrogenation is difficult 
or even impossible. To address this issue, catalysts are commonly added to the metal 
hydride systems, and particularly titanium has been shown to lower hydrogen desorption 
temperatures and enable reversibility |5j, albeit the precise mechanism remains a topic 
of debate (for a discussion, see Refs. [6j [7]). 

Studies of transition metal-catalysed hydrogen storage properties of sodium alanate 
(NaAlH 4 ) went beyond Ti, and also considered Sc, Zr, and other early transition metals 
[H [101 CHI EE2l E3] . Still, a fundamental understanding of how these catalysts, and in 
particular Ti, function has not been achieved yet [HI HHJ [161 E] • It is widely believed 
that a better insight into how a catalyst works could facilitate the rational design of 
novel catalysts for alanates, borohydrides, and amides. 

Carbon nanomaterials have recently been added to the list of materials that could 
potentially act as "true catalysts" in the sense that they retain their structure and are 
not consumed in the reaction, and furthermore would be conceptually much simpler to 
understand in their action [HI HH1 1201 I2H [22j [23]. Besides this function as catalysts 
[211 1251 [261 1271 1281 [291 1301 13T] , carbon nanomaterials have also been considered for their 
role in nano confinement [271 E21 [331 EH EH EH1 E7J [38] to aid in the hydrogen desorption 
process. 

Previously, we presented the results of our study combining experiment and theory 
to demonstrate and understand the catalytic effects of various carbon nanomaterials on 
the hydrogen release and uptake in NaAlH4 |23j. We found that the electron affinity 
of the studied carbon nanomaterials is mainly determined by the surface curvature and 
directly affects the hydrogen sorption mechanism of NaAlH4. One mechanism identified 
by us considers the higher electron affinity causing Na to donate more of its electronic 
charge to the carbon substrate, thus destabilizing the AIH4 complex, which ultimately 
leads to a reduction in the H-removal energy. More recently, the experimental studies 
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were expanded to also cover hydrogen sorption in LiBH 4 , for which it was found again 
that carbon nanomaterials, and in particular fullerenes (Cgo), are excellent in their 
function as a catalyst [39] . 

Hydrogen desorption from LiBH 4 takes place in two steps: 



Step [T] occurs above 400°C and releases 13.8 wt% hydrogen, while step [2] releases 18.4 
wt% hydrogen and takes place at about 950°C which is too high a temperature to be 
considered practical (see Ref. [39] and references therein). 

Here, we present the results from our first-principles analysis of the catalysing 
mechanism that enables Ceo to enhance the hydrogen desorption from LiBH4. In our 
previous work on NaAlH 4 , we had only considered the removal of a single hydrogen atom 
[23] , but in the present study, we went further and explored the stepwise removal of up 
to three hydrogen atoms from LiBH 4 . This approach has led us to identify some new 
features in the Cgo-catalysed desorption reactions of LiBH 4 as compared to NaAlH 4 . 

2. Computational methods 

Total energy calculations were carried out within the framework of the generalized 
gradient approximation (GGA) [30] to density functional theory [3U 132] by using the 
projector augmented wave (PAW) method [43] as implemented in the Vienna Ab initio 
Simulation Package (VASP) [33] ■ PAW potentials with the valence states Is for H, 2s for 
Li, 2s2p for B, and 2s2p for C were employed. The calculations under periodic-boundary 
conditions utilized a cubic supercell box with a side length of 20 A, which prevented 
any unphysical interaction between repeated images and allowed us to consider only 
the T-point for Brillouin zone sampling. Atomic positions were relaxed with respect to 
minimum forces using a conjugate- gradient algorithm. To ensure that we identify indeed 
the ground-state structure, we initiated the relaxation process from several different 
starting configurations which saw the LiBRj unit (for x = 4, 3, 2, and 1) located at 
various potential trapping sites on the fullerene surface (such as on top of a C atom, on 
a C-C bridge, or a hollow site) and in various orientations. 

3. Results and discussion 



Before introducing C 60 , we calculated the energetics for isolated BH 4 _ I units (x = 
0, 1, 2, and 3), both in the charge-neutral state (BH 4 _ rr °) and with one extra electron 
added (BH 4 _ X .~), and for the corresponding isolated LiBH-j. units. From this data, we 
calculated the gain in energy when a hydrogen atom is added and compared the results 
within each of these three systems, i.e., charged and neutral BH 4 _ X and LiBH x . The 
findings from this analysis are summarized in figure [TJ Quite analogous to our earlier 



2LiBH 4 -> 2LiH + 2B + 3H 2 



2LiH -> 2Li + H 2 



(1) 
(2) 



BH 4 7/ 
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published analysis for AlH 4 _ ;c /~ [23J, it is seen that the largest stabilization due to the 
extra electron occurs for the complex containing four hydrogen atoms, i.e., BH 4 . 

We now turn to the main part of our investigation, namely the stepwise 
dehydrogenation in LiBH 4 when supported on C 60 . In the following, we present and 
discuss the results from our density functional calculations for each of the LiBH 4 _a;+C6o 
entities where x takes on the values 0, 1, 2, and 3. 
LiBH 4 

The most stable configuration has Li of the LiBH 4 unit centered above a C-C 
bridge (figure |2^l) . Other configurations which had LiBH 4 located at different sites on 
Ceo were tested, but were found to yield higher energies. For example, a 90°-rotation 
of the LiBH 4 unit (so that both Li and B would have roughly comparable distances 
to the fullerene surface) results in a drastic increase in energy by nearly 1 eV. Also, 
placing the LiBH 4 unit directly above the center of the hexagon yields a state about 
0.25 eV higher in energy. The binding energy of LiBH 4 to C 60 is calculated as 0.50 eV, 
implying that the interaction is rather weak. A Bader analysis of the charge density 
shows that lithium has completely lost its electron. But this charge does not appear to 
have been transferred to the fullerene; rather, LiBH 4 as a whole is charge-neutral, so 
here the electron seems to have been used to stabilize the BH 4 ~ unit, in line with what 
one would expect based on the results from our analysis of the energetics of the charged 
and neutral BH X units above. In addition, partial density of states (PDOS) calculations 
show that the B-H bonds in the BH 4 ~ unit are of covalent nature, while the interaction 
between Li + and BH 4 ~ is of ionic nature. 
LiBH 3 

Upon removal of the first hydrogen atom, we arrive at LiBH3 supported on Ceo- 
Remarkably here, a configuration of the LiBH 3 unit similar to the upright one for LiBH 4 
is not the most stable configuration. Instead, a "sideways" orientation (figure |2b) yields 
the ground state (0.86 eV lower in energy than the "upright" orientation). This enables 
the boron of the BH 3 unit to interact with one carbon atom of the fullerene, thus allowing 
it to preserve an approximate tetrahedral structure in the form of BH3C. One can even 
notice that the carbon atom has been slightly "pulled out" by the boron, distorting the 
hexagon somewhat. Indeed, PDOS calculations show that a bond of covalent nature is 
formed between B and this particular C atom. In a sense, this interaction of BH3 with 
C 60 compensates for the loss of a hydrogen and substantially lowers the energy, as can 
also be seen from the greatly enhanced binding energy of LiBH3 on Ceo (2.09 eV), and 
this plays a crucial role for the low H- removal energy of 2.82 eV in the first step (figure 
|3]). However, it should be noted that the actual energy cost for removing the hydrogen 
atom is higher than that, due to a reaction barrier which was calculated by us with 
the nudged elastic band method [4.5] to be 4.3 eV high. As for the charge state, it is 
seen that lithium retains about 10% of its electron, but the LiBH3 unit as a whole is 
missing —0.7 e, so here a substantial charge transfer to the C 60 has indeed taken place. 
In fact, the carbon atom closest to boron carries an excess —0.37 e. The rest of the 
charge is distributed over the other 59 carbon atoms of the fullerene. Thus, a complex 
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cooperative process is at work here, in which the extra electron that stabilized the BH 4 
unit is no longer needed in BH 3 following the H atom removal, and can now instead 
facilitate the formation of a bond between B and C. 
LiBH 2 

Removal of the second hydrogen atom leads to LiBH2 which just like LiBH3 prefers 
a "sideways" orientation that is found to be more stable than the "upright" orientation 
(figure [2b). The difference in energy between these two possible configurations amounts 
in this case even to 3.5 eV, apparently owing to the now possible stronger interaction 
of boron with two carbon atoms of the fullerene. The binding energy of LiBH 2 to Cgo 
is 4.40 eV, about twice as much as in the case for LiBH3, reflecting that B interacts 
with two C atoms now. This is also seen from the PDOS, showing the formation of 
two covalent bonds between B and the two C atoms. The tetrahedral structure is again 
approximately preserved in the BH2C2 configuration that is formed. Thus the two 
carbon atoms act again as substitutes for the removed hydrogen atoms and are seen 
to be "pulled up", away from the Ceo, towards the boron, without actually breaking 
the Ceo framework. The lowering of the total energy through this form of interaction 
between BH2 and the fullerene is once again responsible for the even further lowered 
H-removal energy of 2.38 eV in the second step (figure [3]). The charge state situation 
is somewhat similar to the previous case of LiBH 3 , in the sense that lithium is charged 
positively with +0.88 e, but the LiBH 2 now carries an even higher positive charge of 
+1.4 e, twice as much as in the previous case. Of course, here, boron is interacting with 
two carbon atoms and indeed it is seen that they both receive about —0.5 e in excess 
charge. 
LiBH 

After the third hydrogen atom has been removed, it is no longer possible for boron 
to easily preserve the preferred tetrahedral arrangement of binding partners as seen in 
the previous three cases. Instead, the BH unit continues to interact with two carbon 
atoms and only to a weaker extent with lithium (figure EH). As a consequence, the 
H-removal energy in the third step is found to be not lowered, but quite oppositely, 
increased with respect to the situation of an isolated LiBH unit where no C 60 is present 
(figure [3]). The binding energy of LiBH to Ceo amounts to 3.66 eV, similar to the value 
for LiBH2, due to both having B effectively interact with two C atoms, but reduced 
perhaps because of the strained non-tetrahedral configuration in which B finds itself. 
The PDOS again shows the formation of two covalent bonds between B and the two 
C atoms. We also tested the possibility of H attaching to Li rather than B, but the 
relaxation results indicate that H does indeed prefer to bind to B in this cluster. 

4. Concluding remarks 

Carbon nanomaterials certainly show great promise as catalysts for hydrogen sorption in 
complex light metal hydrides. In the present theoretical study, we have concentrated on 
the fullerene Ceo acting on the borohydride LiBH4, motivated by definitive experimental 
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evidence for the catalysing effects in this particular combination But as we 

had shown earlier [23J, a similar outcome was seen for the alanate NaAlH 4 , and the 
carbonaceous catalysing agent was found to be not limited to Ceo, but could for example 
as well be high-curvature single-walled carbon nanotubes. Electron affinity certainly 
plays an important role, but not necessarily always the same role. For NaAlH4 on 
various carbon substrates, we had found a charge transfer which presumably destabilized 
the A1H 4 complex [23J. Contrary to that, we did not find any theoretical evidence in 
our calculations here for charge transfer occurring from LiBH-4 to Ceo- Instead, the 
mechanism, which led to a sizeable reduction in the energy required to remove the first 
hydrogen atom in LiBH-4 relies on an increased stabilization of the product state, namely 
by formation of a "substitutional" bond between the B atom of LiBH 3 and a C atom 
of Ceo- That bond is however facilitated by a charge transfer reaction; hence electron 
affinity comes into the equation again. An analogous mechanism, now involving two C 
atoms, is responsible for the even further reduced hydrogen removal energy when going 
from LiBH3 to LiBH2. Those are the main findings of this theoretical investigation. 

Obviously, there are ways in which the present study could be improved. For 
example, our cluster approach only considered a single LiBH4 unit interacting with 
Gqq. Although often, such minimal clusters mirror many of the important qualitative 
properties of the parent solid that they stand for [16], it might still be worthwhile to 
expand the cluster to include more than one LiBH4 unit. Of course, this will significantly 
raise the complexity of the investigation, but could be the only way to identify overlooked 
mechanisms in which a cooperative interaction between several LiBH-4 units on the 
surface of a Ceo leads to a reduction in hydrogen removal energy. Certainly, there are 
many more questions to be answered before we have a complete understanding of the 
catalysing effects of carbon nanomaterials on complex light metal hydrides, and we 
hope that the present work has not merely contributed some new insights to this field, 
but, perhaps more importantly, will also stimulate fresh theoretical and experimental 
investigations. 
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2 3 4 
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Figure 1. This plot shows the energy gain when adding a hydrogen atom to the 
respective clusters of BH„_i (left-side green bars), BH,^_ 1 (central red bars), and 
LiBH„_i (right-side yellow bars). 




Figure 2. Step- wise dehydrogenation of LiBH4 catalysed by Cgo- Shown here are the 
atomic configurations as obtained from ab initio geometry optimizations for fullerene 
plus (a) LiBH4, (b) L1BH3, (c) LiBH2, and (d) LiBH. Li is colored in red, B in green, 
H in white, and C in gray. 
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Figure 3. The hydrogen removal energies are plotted as a function of hydrogen content 
in LiBHa; (from left to right for x = 4, 3, and 2 H atoms), in the absence of fullerene 
(left-side red bars) and when fullerene is present (right-side blue bars). The relative 
change in hydrogen removal energy due to presence of the Ceo catalyst is printed above 
the bars. 



